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ABSTRACT 

The escape fraction, / csc , of ionizing photons from high-redshift galaxies is a key parameter 
to understand cosmic reionization and star formation history. Yet, in spite of many efforts, it 
remains largely uncertain. We propose a novel, semi-empirical approach based on a simulta- 
neous match of the most recently determined Luminosity Functions (LF) of galaxies in the 
redshift range 6 ^ z ^ 10 with reionization models constrained by a large variety of exper- 
imental data. From this procedure we obtain the evolution of the best-fit values of / csc along 
with their 2-a limits. We find that, averaged over the galaxy population, (i) the escape fraction 



0.179 



-0.331 
-0.132 



at z 



8; (ii) at z = 10 



increases from / csc = 0.068lo 047 at z = 6 to / e; 
we can only put a lower limit of / csc > 0.146. Thus, although errors are large, there is an 
indication of a 2.6 times increase of the average escape fraction from z = 6 to z = 8 which 
might partially release the "starving reionization" problem. 
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1 INTRODUCTION 

One of the most crucial issues regarding the evolution of inter- 
galactic medium (IGM) and cosmic reionization is the escape 
fraction, / csc , of ionizing photons from high-redshift galaxies. 
This parameter remains poorly constrained in spite of the many 
theoretical and observational attempts made in past few years. 
The difficulties largely arise from the lack of a full understanding 
of the physics of star formation, radiative transfer and feed- 
back processes, and from uncertainties on the properties of the 
high-z galaxy interstellar medium (ISM); a s a result, derived 
value s of / csc span the large range 0.01 — 1 dFernandez & ShuiH 
1201 lh . Observationa lly, / eac can be rel i ably estimated only at 
redsh if ts z < 3 (Leither er et alj 1 19951: iD ove. Shull, & Ferrara 
200Ct iHeckman et al l l200ll; Iciardi. Bianchi & Ferraral 12003 
Giallongo et al.l 120021: iFernandez-Soto et al.1 120031: llnoue et al 



sive so far, as illustrated by their often conflicting results in terms 
of /esc values and trend with redshift and galaxy mass. 

One key aspect of reionization lies in its close coupling 
with the properties and evolution o f first l uminous sources 
(for r eviews, see |Loeb & Barkanal 200 ll: iBarkana & Loebl 
120011 : IChoudhurv & Ferraral l2006al : IChoudhurvl l2009h . Ob- 
servations of cosmic microwave background (CMB) and 
highest redshift QSOs put very tight constraints on the 
reionization history; these allow to construct self-co n sistent 
models of str u cture f ormation dChoudhurv & Ferraral 1 2005 1: 
Wvithe & Loebl 120051 : iGallerani. Choudhurv & Ferraral I2OO6I: 
Choudhurv & Ferraral l2006bl; iDiikstra Wvithe & HaimarJ 



2005 ; iBergvall et alj 120061 : IShaplev et alj l2006l : [Vanzella et al 
2010). On the other hand, theoretic a l stud ies (Wood & Loeb 
20001: iRazoumov & Sommer-Larser] 120061; iGnedml 1 20011 : 
Gnedin, Kravtsov & Chenl |2008t ISrbinovskv & Wvithd hold: 



20071; ISamui, Srianand & Subramaniad 120071 : llliev et al.1 l200l 
Kulk arni & Choudhurvl 1201 lh . The most favorable model, 



2011; 



Razounio\ & Sommer-Larsenl |2010| : Yaiima. Choi & Nagaminel 
201 it iFernandez & Shulll 1201 it 
have been rather inconclu- 
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Kuhlen & Faucher-Giguere 120 lz 



* E-mail: smitra@hri.res.in 
f E-mail: andrea.ferrara@sns.it 
X E-mail: tirth@ncra.tifr.res. in 



which is consistent with the Thomson scattering optical d epth 
T e i = 0.088 ± 0.015 from WMAP7 data dLarson et alfcoid) and 
the Gunn-Peterson optical d epth evolution from QSO absorption 
line experiments at z > 6 dFan et ai]|2006l) . suggests that reion- 
ization is an extended process over the redshift range 6 < z < 15 
dChoudhurv & Ferraral l2006bl : iMitra, Choudhurv & Ferraral |20 111 
120121) . This model also indicates that reionization feeds back on 
star f ormation by suppressing it in the low-mass haloes a t early 
times dThoul & Weinberdll996l : IChoudhurv & Ferrarj2007l) . 

In parallel, direct observations of galaxies at epochs close 
to the end of reionizati on have made astonishing progresses 
over t he past few years dBouwens & Illingworthl I2OO6I; live et al 



120061 : iBouwens et alj . 2007. 20081; lota et alJl2008tlBouwens et al 
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2009; Henry et al. 2009; Bouwens et al. 2010; Oesch et al. 


201C 


Bouwens et alj 


2010; iMcLureetal. 20ld: besch et alj 


2012 


Bradley et alJl2012b. allowing to derive the galaxy UV Luminos- 



ity Function (LF) up to z m 1( 
iBouwens etaUbOlol ; besch et alj 12012b . and to better constrain 
light production by reionization sources. 

Here we aim at combining data-constrained reionization his- 
tories and the evolution of the LF of early galaxies to get 
an empirical determination of the escape fraction. The study 
also provides relatively tight constr aints also on the evolution 
of the star-forming efficiency e«|Faucher-Giguere et al 2008a; 
iKuhlen & Faucher-Gigueral2012h . Throughout the paper, we as- 
sume a flat Universe with cosmological parameters given by the 
WMAP7 best-fit values: ft m = 0.27, fi A = 1 - O m , Q, b h 2 = 
0.023, and h — 0.71. The parameters defining the linear dark mat- 
ter power spectrum are a$ = 0.81, n s = 0.97, dn s /dlnfc = 
dLarson et alfcOlOn . Unless mentioned, quoted errors are 2a. 



150 



100 




2 DATA-CONSTRAINED REIONIZATION 

We start by summarizing the main features of the semi-analytical 
model use d in this work, which is based on IChoudhurv & Ferraral 
j2005b and lChoudhurv & FerranJ d2006bl) . 

The model follows the ionization and thermal histories of 
neutral, HII and Helll regions simultaneously also accounting for 
IGM inhomogeneities described by a lognormal distribution as in 
iMiralda-Escude. Haehnelt. & Reesl d2000h . Sources of ionizing ra- 
diation are stars and quasars. The stellar sources, all characterized 
by a Salpeter IMF in the mass range M* = 1 — 100A#q , can be di- 
vided into two classes, namely, (i) metal-free (i.e. PopIII) stars; (ii) 
PopII stars with sub-solar metallicities. The transition is based on a 
local critical metallicity criterion. Radiative feedback, suppressing 
star formation in low-mass haloes, is included through a Jeans mass 
prescription based on the evolution of the thermal properties of the 
IGM. 

Given the collapsed fraction / co ii of dark matter haloes, the 
production rate of ionizing photons in the IGM is 

d/coii 



hp h (z) = n b N ic 



dt 



(1) 



where n b is the IGM number density, and JVi on is the number of 
photons entering the IGM per baryon included into stars. The pa- 
rameter iVion can actually be written as a combination of three pa- 
rameters: the star-forming efficiency e„ (fraction of baryons within 
collapsed halos going into stars), / csc , and the specific number of 
photons emitted per baryon in stars, iV 7 , which depends on the stel- 
lar IMF and the corresponding stellar spectrum: 



dv 



dJV„ 



dAf, 



= £*f esc Nj 



(2) 



In our previous work dMitra. Choudhury &FerrariJ[20TlLl2012h . 
we assumed -/Vj on to be an unknown function of z and decompose 
it into its principal components. The Principal Component Analy- 
sis filters out components of the model that are most sensitive to 
the data and thus most accurately constrained. In the following we 
assume a single stellar population (PopII) when computing the ion- 
izing radiation properties; any change in the characteristics of these 
stars over time would be accounted for indirectly by the evolution 
of iVi on . We also include the contribution of quasars at z < 6 as- 
suming that they have negligible effects on IGM at higher redshifts; 
however, they are significant sources of ionizing photons at z < 4. 



Figure 1. Redshift evolution of N- lon obtained from the Principal Compo- 
nent Analysis using WMAP7 data. The solid line corresponds to the model 
described by mean values of the parameters while the shaded region corre- 
sponds to 2-<j limits. 



From the above model, we obtain the redshift evolution of 
iVion by doing a detailed likelihood analysis using three different 
data sets - the photoionization rates Tpi obtained using Lya forest 
Gunn-Peterson optical dept h observations and a lar ge set of hy- 
drodynamical simulations dBolton & Haehneltll2007l) . the redshift 
distribution of Lyman Li mit Systems dN^/dz in 0.36 < z < 6 
dSongaila & Cowiel2010l) and the ang ular power spectra Ci for TT, 
TE and EE modes using WMAP7 dLarson et allhoioh and fore- 
casted PLANCK data. We show the redshift evolution of A'ion(z) 
obtained from our Principal Component Analysis using WMAP7 
data in Fig.[T] The solid line corresponds to the model described by 
mean values of the parameters which we obtained by performing 
a Monte-Carlo Markov Chain (MCMC) analysis over the parame- 
ter space of our model, while the shaded region corresponds to its 
2-cr limits. We concluded that it is not possible to match available 
reionization data with a constant Nj on ove r the whole redshift range 
dMitra. Choudhury & Ferrarj201 11120121) . Rather, it must increase 
at z > 6 from its constant value at lower redshifts. This is a sig- 
nature of either a varying IMF and/or evolution in the star-forming 
efficiency and/or photon escape fraction of galaxies, as eq.[2]clearly 
shows. 

At this point, it is worth pointing out some of the caveats in our 
formalism based on PCA. The MCMC analysis for this work wa s 
done using 2-8 PCA modes dMitra. C houdhury & Ferrar3l2012h . 
Truncating a PCA expansion can reduce the variance in the esti- 
mation of the reionization history, but also introduces a bias to- 
wards the fiducial history. In order to account for this, we used 
the Akaike information criteria (AIC) to reduce the intrinsic bias 
present in any particular choice of fiducial model. We found that at 
z ^ 6, the strong Lyman-a forest constraints essentially fix Mon, 
so that the efforts to reconstruct the reionization history are very 
promising at these redshifts. On the other hand, it is much more 
difficult to recover ]Vi on or the other various quantities related to 
reionization history at z > 6 in a truly model-independent manner 
as the re exists a considerable amount of bias at this high redshift 
end dMitra, Choudhury & Ferrardl201ll . l2012h . It is possible that 
for these redshifts, the statistical uncertainties may have be en hid- 
den as systematic uncertainties dHuterer & Starkmanll2003h . How- 
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ever, with more data it would be possible to apply this technique 
in a regime where the variance in iVi on is small enough to produce 
a useful constraint on the reionization history without the need to 
truncate the PCA modes so severely and so without introducing 
any bias. This technique will become more applicable as more data 
becomes available for z > 6 region. 



3 LUMINOSITY FUNCTION EVOLUTION 



The 
LF 



effect 
was 



of reionization 
studied using 



on the high redshift galaxy 
the semi-analytical mod- 



els by Samui. Srianand & Subramanianl ( 120071) and 



iKulkarni & Choudhurvl 1 201 1 ). In this work, we follow their 



method to study the evolution of LF for our model. 

The LF is derived as follows. We compute the luminosity at 
1500 A of a galaxy having the halo mass M and age At using 



i.r.M„(.\/.Arl = r ( ) Ml ir (Afi 



(3) 



Here the age of the galaxy formed at z' and observed at z is At = 
t z ~ t z i , iisoo(At) is a template specific luminosity at 1500 A for 
the stellar population of age At. As we restrict to a single stellar 
population, i.e. PopII stars, e„ indicates the star forming efficiency 
of PopII stars throughout the paper. 

To compute £1500, we use stellar population models of 
iBruzual & Chariot! d2003h for PopII stars. The UV luminosity de- 
pends on galaxy properties including th e IMF, star f o rmatio n rate 
(SFR), stellar metallicity (Z ) and age. iDaval et all J2009al) and 
Daval. Ferrara & Sard ( feOlOh have shown that the metallicity cor- 
relates with stellar mass, and the best fit mass-metallicity relation 
they find is 



Z/Z @ = (0.25 - 0.05A«) log 10 (M») - (2.0 - 0.3Az) 



(4) 



where Az = (z— 5.7) and M„ is the total stellar mass of the galaxy. 
We take all the available stellar population models in the metallicity 
range Z — 0.0001 — 0.05 for PopII stars and interpolate them to 
compute £1500 following the mass-metallicity relation given by the 
above relation for our model galaxies. 

The lu minosity can be converted to a standard absolute AB 
magnitude ( Oke & Gunn 1983; Sam ui. Srianand & Subram anian 
2007: Kulkarni & Cho udhurv 20 lj) using 



Mab = -2.51o gl 



erg s _1 Hz _1 



51.60 



(5) 



The luminosity function Q(Mab, z) at any redshift z is then given 

by 

dn An dLi 50 o 



AB,Z 



dM A i 



dZ/1500 dM A B 



where 

dn 
dLi5oc 



Az' 



AM 



(L1500, At) 



AMAz 



7 (M,z') 



(6) 



(7) 



is the comoving number of objects at redshift z with observed lumi- 
nosity within [Li5oo,ii50o + dLisoo]- The quantity d 2 n/dMdz' 
gi ves the formation rate o f halo es of mass M, which we obtain as 
in Ichoudhurv & Ferraral d2007l) . Note that, we can vary the star- 
forming efficiency e* in eq. [3] as a free parameter and obtain its 
best-fit value by comparing the high-redshift LFs computed using 
the above equations with observations. 
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Figure 2. Luminosity function from our model for best-fit e* (black curve) 
and its 2-cr limits (shaded region) at z = 6, 7, 8 and 10. Data points with 
2-cr errors are from Bouwens & Illingworth (2006) (z = 6), Bouwens et al. 
(2010) (z = 7, 8) and Oesch et al. (2012) (z = 10). For 2 = 10, we show 
the luminosity function from our model for the maximum value of e* for 
which the LF curve does not exceed the experimental upper limits. 



3.1 Constraining the escape fraction 

Our strategy to constrain / osc exploits the combination between the 
previously derived (Sec. [2} evolution of Ai on , and the constraints 
on e» that can be derived from matching LFs at different redshifts. 
Once the (Salpeter) IMF of the (PopII) stars is fixed, is also 
fixed and equal to w 3200; from eq.[2]we then get the value of / csc 
as follows: 



A io 



(8) 



As the uncertainties on [Mon/iV 7 ] and e, are independent, the 
fraction al uncertainty in / esc can be obtained from the quadrature 
method jTavloifT9 97)). i.e. 



& fcsc 
fesc 



S [A ion /A 7 



+ 



(9) 



In this work, we are interested in the z ^ 6 evolution of the escape 
fraction. In principle, our approach can also be used for the lower 
redshift range 3 ^ z ^ 5, provided that a detailed treatment of dust 
extinction is added to our model. The underlying assumption in the 
present work is that dust effects on the escape fraction can be safely 
neglected at early times. 



4 RESULTS 



The observationally 



deter mined 

iBouwens & Illingwortr] (2006) for 



LFs 



fcOlCh for z 



re taken from 
6, iBouwens et al.l 
7, 8 and lOeschet"ai] ( l2012l) for 2 = 10. Figured] 



shows the globally averaged LFs calculated using our model for 
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Redshift 


Best-fit e* 


2-cr limits 


Best-fit /esc 


2-cr limits 


2 = 6 


0.0365 


[0.0253,0.0481] 


0.0684 


[0.0210,0.1221] 


2 = 7 


0.0385 


[0.0193, 0.0576] 


0.1607 


[0.0319, 0.4451] 


2 = 8 


0.0523 


[0.0129, 0.0822] 


0.1794 


[0.0466,0.5098] 


2 = 10 


< 0.0841 




> 0.1456 





Table 1. Best-fit values and 2-cr limits of e» and the derived parameter / OS c for the reionization model obtained from the LF calculation at different redshifts. 
At 2 = 10, we get only an upper limit of e„ and a corresponding lower limit of / eS c. 
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Figure 3. Redshift evolution of the escape fraction / csc with 2-cr errors. The 
2 = 10 data point shows the lower limit of / csc . The solid line corresponds 
to its best-fit value while the shaded region corresponds to 2-cr limits. 



z — 6, 7, 8, 10 compared to the observational data points. The 
z = 10 data are obtained from the detection of a single galaxy 
candidate by lOesch et al.l (2012); hence, we only show results for 
the maximum value of e» for which the LF curve does not exceed 
the experimental upper limits. 

Our model reproduces the observed LFs reasonably well, es- 
pecially at lower redshifts. From such a match we find that the best- 
fit value of the star-formation efficiency e* nominally increases 
from 3.6% at z — 6 to 5.2% at z = 8. Such a small variation is 
statistically consistent with a constant value of e„, i.e. no evolution. 

The corresponding values of / csc calculated using eq.[8]and|9] 
are plotted in Fig. [3] along with the 2-cr confidence limits (shaded 
region). The numerical values for e„ and / csc are also reported in 
Table Q] for different redshifts (z = 6,7,8). The escape fraction 



shows a moderately increasing trend from f C! 



0.068 



+0.054 



at 



6 to / c . 



0.179 



+0.331 



at z 



at z = 10 we can only 



put a lower limit of / csc > 0.146, corresponding to the maximum 
allowed value of e* = 0.0841. 

The reported 2-cr errors are however relatively large and we 
cannot exclude a non-evolving galaxy-averaged trend for / eS c • The 
uncertainties become larger with redshift as a consequence of the 
fact that the larger LF errors at higher redshifts. 



5 CONCLUSIONS 

We have used a semi-a nalytical model, based on 

IChoudhurv & Ferrar A b005h and lChoudhurv & Ferrar <J fe006bl) to 
compare the UV luminosity functions at different epochs predicted 



from our model with the observed LF to constrain the parameters 
related to star formation history in the redshift range 6 ^ z ^ 10. 
In particular, by varying the star formation efficiency as a free 
parameter, we have constrained one of the most unknown param- 
eters of reionization models, the escape fraction / csc of ionizing 
photons from high-redshift galaxies. The main findings of our 
work are that, averaged over the galaxy population, (i) the escape 
fraction shows a moderate increase from / csc = 0.068tn'rS* at 



6 to / csc = 0.179^ 



at z 



-0.047 

(ii) at z = 10 we can 



only put a lower limit of / csc > 0.146. Thus, although errors are 
large, there is an indication of a 2.6 times increase of the average 
escape fraction from z — 6 to z = 8 which might partially release 
the "starving reionization" problem. At the same time, the best-fit 
value of the star formation efficiency e* nominally increases 
from 3.6% at z — 6 to 5.2% at z = 8. Such a small variation is 
statistically consistent with a constant value of e», i.e. no evolution. 

Parallel to our more phenomenological approach, in the past 
few years many numerical and analytical studies have attempted to 
constrain / csc reaching often contradictory conclusions, likely due 
to uncertainties on star formation history, feedba ck, radiation trans 
fer an d the geometry of the ISM distributio n ( Fernandez & Shul 
2011 ). Increasing iRazoumov & Sommer-Larsenl |2006| , 1201 



Haard t & Madaul 20111), decreasing dWood & Loebl|2000l) or un- 
evolving jGnediri 12008 : lYai ima. Choi & Nag amine 2011) trends 



have been suggested as a function of redshift. 

A s trong redshift evolution of the escap e fraction was recently 
found by Kuhlen & F aucher-Giguerej J20 1 2). They show that, mod- 
els in which star formation is strongly suppressed in low-mass 
haloes, can simultaneously satisfy reionization and lower redshift 
Lyman-a forest constraints only if the escape fraction of ionizing 
radiation increases from ~ 4% at z — 4 to ~ 1 at higher red- 
shifts. Although broadly in agreement with their conclusions, our 
results show instead that reionization and LF data can be satisfied 
simultaneously if / csc grows from ~ 7% at 2 = 6 to ~ 18% at 
z — 8, but without requiring an escape fraction of order of unity at 
these redshifts. We believe that t his discrepancy can be understoo d 
as due to the fact that unlike Kuhlen & Faucher- Giguerd J2012h . 
we are fitting the/;;// CMB spectrum rather than the single value of 
T e v, the latter choice can be thought as a simplification of CMB po- 
larization observations. In addition, we have used a PC A analysis 
to optimize model parameters to reionization data, yielding a more 
robust statistical analysis. 

Although here we have only considered the evolution of z ^ 6 
luminosity functions, our approach can also be applied to model 
the LFs at 3 ^ z ^ 5. As hydrogen reionization mostly occurs at 
z > 6, the LFs in this lower redshift range are very unlike to be sen- 
sitive to the details of reionization history. Also, dust extinction at 
2 < 6 c an decrease / esc by absorbing t he ionizing photons at these 
epochs jYaiima. Choi & Nagamine 201 1). As a caveat we mention 
that the present results can be responsive to changes in some c os- 
mological parameters, mainly as and n s {Pandolfi et al. 2 0"llT) . A 
larger erg or n s may lead to an increase in the number of collapsed 
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haloes at all redshifts. In principle then, one should include these 
two quantities in the analysis as additional free parameters. Also, it 
could be interesting to evaluate the effects of PopIII stars and other 
feedback processes in our LF calculation. We hope to revisit some 
of these topics in more detail in future work. 
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